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Abstract 

The mean-square displacement (MSD) was measured by neutron scattering at various tempera- 
tures and pressures for a number of molecular glass-forming liquids. The MSD is invariant along 
the glass-transition line at the pressure studied, thus establishing an "intrinsic" Lindemann crite- 
rion for any given liquid. A one-to-one connection between the MSD's temperature dependence 
and the liquid's fragility is found when the MSD is evaluated on a time scale of ~4 nanoseconds, 
but does not hold when the MSD is evaluated at shorter times. The findings are discussed in 
terms of the elastic model and the role of relaxations, and the correlations between slow and fast 
dynamics are addressed. 
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I. INTRODUCTION 



A major challenge to current condensed-matter physics is to explain properties of su- 
percooled liquids approaching the glass transition. In particular, there is still no consen- 
sus about what causes the dramatic " super- Arrhenius" increase of the liquid's viscosity. 
The viscosity will for so-called fragile liquids close to the glass transition increase with 
about one decade (in some cases even more) for a temperature decrease of just 1%. The 
most popular model is the Adam-Gibbs entropy modeli which relates the slowing down 
to an underlying phase transition; this model involves the notion of a growing length 
scale as the liquid is cooled towards the glass transition. The experimental search for 
growing length scales uses many different observables, it has been going on for several 
decades^ and it is an ongoing active field of research^"-, but there is still no clear con- 
clusion. Another approach is advocated by the so-called elastic models, the main ideas 
of which date back to a paper by Eyring and co-workers from 19432 and Nemilov's fur- 
ther work in 1968^. These models relate the activation energy to the viscous liquid's 
short-time elastic properties^^i^. Thus, apparently paradoxical, properties on the pico- or 
nanosecond time scale could determine the slow molecular relaxations taking place over 
minutes or hours in the liquid close to the conventional glass-transition temperature. 
This paper investigates this intriguing prediction by studying several liquids, some un- 
der varying pressure, by means of quasi-elastic neutron scattering experiments providing 
the most direct measurement available of the mean-square displacement on several short 
time scales. 

The idea of a cormection between the dynamics on time scales differing by ten or more 
orders of magnitude has also been put forward in other contexts. If correct, it emphasizes 
the fact that the glass-transition phenomenon involves an exceedingly large dynamical 
range — and that a full understanding of the glass transition must encompass both fast 
and slow dynamics. Before proceeding to describe the experiment, we give examples of 
the ideas and results which point in this direction. 

In 1992 it was observed by Buchenau and Zorr>ii that there is a relation between the 
temperature dependence of the structural relaxation time and the temperature depen- 
dence of the mean-square displacement (MSD) observed on the nano-pico second time 
scale in selenium as determined by neutron scattering experiments. Since then, other 
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groups found similar results, showing qualitatively that the larger the fragility is, the 
stronger is the temperature dependence of the MSD^~— . Recent theoretical work relates 
the mean-square displacement and the vibrational entropy and relates both quantities to 
the slow dynamicsi^. In 1987 Hall and Wolynes^Si theoretically discussed how the mean- 
square vibrational displacement controls the relaxation time according to the expression 
T a exp (Const./ (u^))- Their approach was later developed into the random first-order 
transition theory (RFOT) of the glass-transitiork^i*^^ where a variational density profile 
built of Gaussian vibrational displacements around aperiodic atomic positions is opti- 
mized for free-energy minimization. Thus the vibrational short-time displacement is the 
crucial quantity for the RFOT, which was later developed into a full-fledged theory lead- 
ing to a generalization of the Adam-Gibbs prediction for the relaxation time in terms of 
the configurational entropy-. 

In more recent works22<2i Novikov and Sokolov demonstrated a surprising connection 
between the "fragility" — a measure of how fast the liquid's viscosity (or relaxation time) 
increases as temperature decreases and enters the glassy state — and elastic properties 
of the glass: The more fragile the liquid is, the higher is the ratio between the bulk and 
shear moduli of the resulting glass. However, a study of a larger set of liquids shows that 
the correlation does not hold on general^. Novikov and Sokolov discussed a possible 
explanation of their correlation in terms of elastic models like the shoving model^. This 
does not seem to be a correct connection because recent experimental work by Nelson 
and coworkers supports the shoving model while it lends no support to Novikov's and 
Sokolov's correlation^^. Sokolov earlier introduced a parameter derived from the mea- 
sured dynamic structure factor, which relates the strength of the quasielastic scattering 
intensity at Tg normalized to the intensity of the boson peak and the fragility^^. This cor- 
relation has also been questioned^^ even if it holds on a very qualitative level (intense 
boson peaks are seen in strong systems like oxide glasses whereas the boson peak in van 
der Waals systems is less pronounced). It also appears that this correlation could be more 
related to isochoric than isobaric fragility^^. Along similar lines Scopigno and coworkers 
suggested a correlation between the temperature dependence of the non-ergodicity fac- 
tor measured in the glass and the fragility of the liquid^S*^. The main point to be noted 
here is that the two Novikov-Sokolov correlations as well as the Scopigno correlation 
provide an intriguing connection between short- and long-time properties of the liquid 
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(since the liquid's short-time mechanical properties are those of the glass corresponding 
to the liquid structure at the temperature in question). 

Widmer and Harrowell^^ proposed studying the Maxwell-Boltzmann ensemble aver- 
aged mean-squared displacement of a particle for any given initial configuration (the 
"iso-configurational ensemble"), terming this quantity the dynamic propensity of the 
particle in question. This property's distribution reflects the dynamic heterogeneity of 
the liquid. Thus, once again, a connection is established between the long-time dynamic 
properties and MSD on time scales much shorter than the relaxation time. 

Leporini and collaborators^ argued from simulations that there is a universal corre- 
lation between the structural relaxation time and the "rattling amplitude" from high- to 
low-viscosity states. According to this picture the glass softens when the rattling ampli- 
tude exceeds a critical value. This implies a "universal" Lindemann criterion for the glass 
transition, i.e., that the glass transition takes place when the MSD reaches a certain value 
(see e.g. Ref. HQ. 



At first sight it appears very surprising that there could be any relation between the 
alpha-relaxation process — taking place on the second or hour time scale — and mean- 
square displacements taking place on the nano/pico second time scale. It should be re- 
called, however, that whereas the alpha relaxation is very slow, the barrier transitions 
themselves are fast. This fact is the starting point for the elastic models. In these models 
the relation between mean-square displacement and fragility comes very natural, and the 
stiffness of the material or, equivalently, the steepness of the energy minima, determines 
the activation energy of the alpha process^^^^^. 

In terms of the vibrational MSD {u^),m the simplest version where the instantaneous 
bulk and shear moduli are proportional in their temperature variation, all elastic models 
imply for the activation energy AE(T) a Ta^ / {iP-){T) at atmospheric pressure, where a 
is the average intermolecular distancei2*1^20*^. This implies 



T(r) = TO exp . (1) 

where C is a constant. This result relates a larger MSD to a shorter alpha relaxation time 
and it implies that the temperature dependence of the relaxation time is governed by the 
temperature dependence of the MSD. Hence the elastic model predicts that the change of 
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the MSD just above Tg is more dramatic the more fragile the liquid is. This agrees with 
the generally observed trend, although questions remain about at which time scale (u^) 
should be considered. The glass-transition Lindemann criterion states that the (u^) /a^, 
i.e., the relative vibrational amplitude of the atoms, at the glass transition should reach a 
certain universal value allowing diffusion on long time and length scalesi2*2iiM42. Recall 
that the Lindemann criterion is the rule that (u^) / ^1% when any crystal melts. If the 
glass transition is also characterized by such a universal number, there would be an ap- 
pealing analogue between crystal and glass "melting" - although the latter phenomenon 
is known to be cooling rate dependent. 

In this paper we present a quantitative test of Eq. ^ based on MSD data in the nano- 
and pico-second time scale obtained by neutron backscattertng and time-of-flight tech- 
niques on several molecular liquids, covering fragilities ranging from 49 to 161. Three 
liquids were also studied under varying pressure. This allows one to examine different 
glass transitions of the same liquid, i.e., without changing the intermolecular interactions. 

II. EXPERIMENTAL 

The experiments were carried out on the back-scattering instruments INIO and INI 6 at 
the ILL. This spectrometers use the (1 1 1) reflection of Si-single crystals with a Bragg an- 
gle of 90° at the monochromator and analyzers to reach an energy resolution of FWHM=1 

o 

^eV (corresponding to a time scale of ^4 ns). The wavelength of the neutrons was 6.27 A. 
The wavevector, Q, range covered was 0.2A~^ to 1.9 A~^. The experiments were per- 
formed isobarically in cooling with a rate of ^ 0.5 K/ min. Pressure was applied using 
a clamp pressure cell mounted on the bottom of an insert to the cryostat. Sample trans- 
mission was 88% for the high-pressure measurements and 95% for the measurements 
at atmospheric pressure in standard aluminum cells. The liquids studied are: glycerol, 
cumene, dibuthylthalate (DBF), m-toluidtne, sorbitol, triphenylphosphite (TPP) and dec- 
ahydroisoquinoline (DHIQ). The first three were also studied at elevated pressure (300 
MPa or 500 MPa); the pressure dependence of the glass-transition temperature was ob- 
tained from calorimetric experiments or extracted from dielectric spectroscopy under 
pressure. For glycerol a different pressure cell, dedicated for studying liquids under 
hydrostatic pressure, was used. This ceU is built out of Niobium resulting in a low back- 
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ground and the pressure can be adjusted via a capillary from the outside of the cryostat, 
but the maximum pressure is limited to 300 MPa^. The raw data correction was per- 
formed using the standard ILL software Sqwel which converts the measured data to the 
scattering law S{Q,co) for the sample. The MSD is calculated from the measured elastic 
intensities by adopting the Gaussian approximation, ln(I) = A — ^ g" ^ with A being a 
constant. The measured intensity is normalized to the intensity measured at low temper- 
ature, T = 4 K, which means that the zero point motion is removed. We find that the 
dependence is obeyed in the temperature range K to 1.2 Tg^^. Even so, one should 
still be aware that the measured MSD can contain local relaxations which are unrelated 
to the vibrations and independent of the structural relaxation. 

Supplementary experiments on DHIQ and DBF were carried out on the back- 
scattering instrument IN13 at atmospheric pressure. The energy resolution on IN13 is 
almost ten times wider, FWHM=8 }ieV, meaning that the MSD we access in the measure- 
ment is probed on an almost 10 times faster time scale (^0.5 ns). 

III. THE LINDEMANN CRITERION 

Our measurements give the molecular MSD at the nano and the picosecond time scale. 
How can the elastic model prediction be tested? Comparing data for the same liquid at 
differing pressures, avoids making assumptions about the constant C of Eq. [H If it turns 
out that C is common to all liquids, a universal (i.e., genuine) glass-transition Ltndemann 
criterion is implied. 

In Fig. [T]we present data for three liquids of different chemical nature, intermolecular 
interactions and fragility (glycerol, cumene, and dibutylphthalate). They are studied at 
ambient as well as at high pressures (300 MFa for glycerol, 500 MFa for the two other 
liquids). The left part of each figure gives the mean-square displacement as a function of 
temperature at the two pressures where the dashed line marks the glass-transition tem- 
perature. The right part gives the data scaled as implied by Eq. Thus the temperature 
is scaled by Tg and the MSD by a^, where we assume that the intermolecular distance 
scales with oc p"^^^. The density, p is evaluated at {Tg{P),P) from known equations 
of stated. For all three liquids there is data collapse, showing that a Ltndemann type 
criterion is fulfilled. 
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FIG. 1: (Color online) The MSD of glycerol, cumene and DBP at atmospheric pressure and at 
500 MPa (300 MPa for glycerol). The left hand side of the figure shows {u^) and temperature on 
an absolute scale. The dashed lines indicate Tg, the dash-dotted lines Tg at high pressure. The 
temperature scale in the right hand side of the figure is scaled by the pressure dependent Tg and 
the y-axis is scaled with cx p^^/-^ evaluated at {Tg{P), P). 

The scaling of the temperature axis is by far the most important for this data collapse 
in the pressure range studied. The estimated increase in density is less than 10%. This 
gives a decrease of fl^ by approximately 5%. This difference is almost indistinguishable 
in figure [l] due to the scatter of the data. A scaling with the pressure-dependent glass 
transition temperature was earlier shown by two of us for a polymer sample in ref. |45 
The earlier scaling is also done at relatively low pressures where the density does not 
change dramatically. It thus appears that for the systems studied so far, the MSD is con- 
stant along the glass-transition line in a P-T-diagram, suggesting a Lindemann criterion. 
To verify whether the MSD is constant or, as suggested by the elastic model, if the nor- 
malized MSD, {u^){T)/a^, is constant studies are needed in a larger pressure range. A 
constant value of {u^){T)/a^ with oc p~'^^^ along the glass-transition line, CTa{P),P), 



37 



. Alter- 



is also consistent with the existence of isomorphs as it is pointed out in Ref. 
natively, the change of the intermolecular distance a could be found based on a more 
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local measure extracted from the static structure factor, S(Q). However, measurements 
as a function of temperature and pressure show that the Q-dependence of the peak max- 
imum follows the same behavior—. 



FIG. 2: (Color online) The temperature dependence of {ip-) scaled to = v^^^ (see text for details) 
for 5 different liquids; glycerol (circles), DBF (diamonds), m-toluidine (squares), cumene (trian- 
gles), DHIQ (crosses). The temperature is scaled to Tg. The inset: The same data as in the main 
figure, here shown with the absolute value of (w^). 

The next step is to investigate whether the constant C in Eq.HJis common to all liquids, 
as required by a universal glass-transition Lindemann criteriork^*^. This is investigated 
in Fig. |2]by plotting {u^){T) / as function of T/Tg for a selection of liquids at ambient 
pressure with quite similar Tg's. If the constant C were universal, / (u^) should be the 
same for all liquids at Tg. The figure shows that this is not the case since the number 
{u^){Tg)/a^ varies a factor of 3 going from glycerol to m-toluidine. It should be noted, 
though, that the temperature dependence of the MSD of m-toluidine has a strong in- 
crease far below Tg. This type of behavior has earlier been seen in other systems and 
is associated with the methyl-group rotation-^. Such type of local motion is probably 
independent of the glass-transition temperature and an irrelevant contribution to the ap- 
parent MSD with respect to the Lindemann criterion. Even for the four other liquids 
there is a factor 2 in variation when comparing (u^) (Tg) / . Based on these data and as- 
suming that for DBF, DHIQ and cumene the major contributions to the MSD are arising 
from displacements which are relevant for the structural relaxation near Tg, we cannot 
confirm the existence of a universal Lindemann criterion as predicted by Leporini and 
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collaborators'^^, at least not on the time scale we explore. The question is also discussed 
in a recent publication by the same group^. 

IV. TEMPERATURE DEPENDENCE OF THE MSD ABOVE Tg 

We now apply a different scaling for the MSD by normalising to the MSD at Tg. Figure 
|3] shows {u^){T)/ {u^)t„ as a function of T /Tg- The {u^) value of the very fragile liquid 
DHIQ at the nanosecond rises most, the {u^) of glycerol least, dramatically; the three 
remaining liquids, which all have similar intermediate fragilities, fall in between. The 
systems studied hence confirm the general trend that more fragile liquids have more 
temperature dependent amplitude of the short time MSD above Tg than do less fragile^^. 
The elastic models make a quantitative prediction regarding the relation between the 
temperature dependence of {u^) and that of the alpha relaxation time. Thus the elastic 
model leading to Eq. dT) is based on 
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FIG. 3: (Color online) The measured (w^) scaled to {u^)t^ for 5 different liquids. The temperature is 
scaled to Tg. glycerol (circles), DBF (diamonds), m-toluidine (squares), cumene (triangles), DHIQ 
(crosses). The lines above Tg illustrate the fitted slopes used in Fig. ID 
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, it follows 



that the elastic models predict (where the weak temperature dependence of a at constant 
pressure is ignored) 



Using the general relation— between the conventional fragility index and Jp it follows 
that the model predicts a proportionality between Angell's (isobaric) fragility wp and the 
relative change of (w^) with relative change in temperature: 



where tq = 10~ s is the microscopic time and T^=100 s is the relaxation time at the glass- 
transition temperature (where fragility is evaluated). Hence the elastic model predicts a 
correspondence between the slope seen in Fig. |2] at and the fragility found from the 
temperature dependence of the alpha relaxation time. 

Figure m tests this relation using fragilities and T^'s taken from literature (see table HVl 
for values and references). The value of ^^^^j |p(T = T^) is in all cases calculated in the 
temperature range from Tg to ^1.1 T^, corresponding to the range where the fragility is 
determined. The data taken on the nanosecond time scale all lie close to the line. This 
result is rather convincing, especially because Eq. (ID not only predicts that there is a 
proportionality between mp and '^d"n'r^ \ ~ '^s)' value of the proportionality 

constant as well. 

V. THE ROLE OF RELAXATIONS AND ANHARMONICITY 

Although Fig. m shows an overall agreement with the elastic model prediction, when 
using the MSD on the nano second time scale, a number of issues remain to be consid- 
ered. Not only, of course, is a more extensive study of different liquids needed, there are 
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FIG. 4: (Color online) The value of 



dlnT 



as a function of the isobaric fragility. Full symbols 



(blue) are data that refer to the nanosecond time scale (~4 ns), open symbols (red) are data ob- 
tained on IN13 on a ten times faster time scale (~0.4 ns). The blue full line shows the expression 
Eq. (H) which follows from assuming AE(T) cx T/ (u-^) as is predicted from the elastic model (using 
Tg = 100 s and tq = 10^^^ s.). Values and references are given in tableHV] The points calculated at 
~4 ns fall close to the predicted line, while the points calculated at ~0.4 ns where the (w^) is less 
influenced by relaxations fall much below the line (the red dashed curve is a guide to the eye). 
When points lie below the blue full line, the temperature dependence of (m^) underestimates the 
temperature dependence of the activation energy. It should be noted that both x and y variables 
have considerable uncertainty as they are arrived at as numerical derivatives of data. The error 
bars on the literature data are our estimate based on experience and on the difference in reported 
values from different sources. 

also other more fundamental issues. One problem is that in elastic models it is usually 
assumed that the measured (u^) is purely vibrational, i.e., that no relaxational motion 
contributes to (u^) around Tg. It is not likely that this assumption is generally correct, 
how^ever. Thus we know from time-of-flight spectra that DHIQ has a strong quasi-elastic 
scattering already at Tg^^. Time-of-flight measurements have a broader resolution func- 
tion, and consequently shorter time scale, so this quasi-elastic scattering corresponds to 
relaxation at even shorter times than the MSD probed by backscattering. The alpha re- 
laxation itself also enters the experimental window at some time, possibly already when 
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TABLE I: Values and references for the points shown in figure 4. The asterisk indicates that the 
value has been calculated from data in the paper. 

Ta ^ Ijis if the relaxation function is very stretched. This happens intrinsically faster 
for fragile liquids than for strong liquids (for which, also, the relaxation functions are 
generally less stretched). 

When considering relaxation it also appears that the finding of Fig. |4]is consistent with 
another phenomenological feature observed in the dynamic structure factor as measured 
from inelastic scattering. Namely the observation that the relative strength of the boson 
peak compared to the fast relaxation, measured at Tg, is related to the isobaric fragility 
of the glass former: a parameter, R, is defined as the quasielastic intensity divided by 
the boson peak intensity, and proposed to increase with increasing isobaric fragility^^*^. 
This model-independent assertion made by comparing the behavior of different glass 
formers is controversial-^, but appears more robust than many other correlation s ^^ '^^'^^ . 
Our findings are clearly consistent with this observation, showing the importance of fast 
processes on a time scale of a few nanoseconds even close to Tg, where R is determined 
for the most fragile liquids. This observation at the same time suggests that the boson 
peak intensity itself is not the relevant quantity for the correlation, but that it probably is 
the larger intensity of fast relaxation in fragile liquids that yields the correlation. 

To investigate the role of relaxations further we have performed supplementary mea- 
surements of the MSD of DHIQ and DBF ( for glycerol we used literature data^^ referring 
to the same timescale) using IN13 which has a broader than INIO resolution and there- 
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FIG. 5: (Color online) The MSD of DHIQ measured by IN13 (on a time scale of ~0.4 nsec) and 
INIO (on a time scale of ~4 nsec). The dashed lines indicate Tg. The lines above Tg illustrate the 
fitted slopes used in Fig. ID The temperature dependence of the MSD is time scale independent 
below Tg while it becomes strongly time scale dependent above Tg. This indicates that {u^) below 
Tg is dominated by vibrations while relaxations play a role on the nanosecond time scale in the 
liquid above Tg (even close to Tg where the alpha relaxation time is still on the order of seconds). 

for accesses {u^) on a time scale which is approximately 10 times shorter. In Fig. |5]we 
compare the mean-square displacement of DHIQ found on the two different instruments, 
two distinct timescales. The measured {u^) follow each other below Tg, which strongly 
indicates that we probe genuine vibrations in this regime and therefore that the finding 
of the Lindemann criterion is related to the vibrations, as predicted by the elastic models. 
Above Tg, on the other hand, we see a separation of the two curves. It is evident that 
the temperature dependence of the MSD on the nanosecond timescale probed by INIO 
is much more pronounced than the temperature dependence on the shorter time scale 
probed by IN13. This dependence on the timescale indicates that we are not probing 
the purely vibrational MSD on the nanosecond time scale, but rather a combination of 
vibration and fast relaxations. 

To illuminate how the relaxations affect the value of d\r\{y?-) / d\r\T we adopt a 
simple "jump-diffusion" type modelling^-: If on the time scale set by the experiment 
some molecules vibrate whereas others jump once or more, the MSD separates into 
two contributions: {ip-) = {u^)v[h + {u^)\ump- For the log-log derivatives one finds 
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dln{u^) / d\nT = Adln{u^)yii,/dlnT + B d\n{u^)^^^p/d\nT where A = {u^)vih/ W^) 
and B = (H^)jump/ (w^) give the relative weights of the two contributions (A + B = 1). 
The jump contribution is most likely strongly temperature dependent. Thus any cor- 
rection for this in order to get the pure elastic contribution to d\n{u^) / dlnT pushes the 
points in Fig. |4] downwards, i.e., further away from the line. This is exactly what we 
see for the d ln{u^) / d\nT obtained from the IN13 data on DHIQ, which is also shown in 
Fig. m Similarly we see that the d\n{u^) / dlnT calculated from glycerol data taken on 
IN13 reported by Wuttke^ lie below the line. However, the difference between the two 
timescales is much less dominant for glycerol than for the very fragile DHIQ, indicating 
that the relaxation are more dominant in the latter. 

Based on the above considerations, we conclude that using the vibrational part of the 
MSD as deviced by the elastic models (Eq. dT)) underestimates the temperature depen- 
dence of the activation energy. The reason for this could be that Eq. [1] is based on a 
simplified reasoning that basically ignores anharmonicities^. There are two nontrivial 
assumptions going into this reasoning: a) The harmonic approximation, according to 
which the curvature at the minimum is inversely proportional to (w^); b) The energy bar- 
rier being proportional to the curvature at the minimum, as it would be if the potential 
were parabolic and scale accordingly. The first approximation applies to a good approx- 
imation at sufficiently low temperatures and may well apply to highly viscous liquids 
because these have fairly large energy barriers. The second approximation implies that 
if the barrier goes to zero, so does the curvature. However, this is not necessarily the 
case. Thus the simple elastic model assumption that the barrier scales with curvature 
may break down. In summary, for a given temperature dependence of (u^) a more real- 
istic model might well predict larger d In AE/d In T than predicted by the elastic models. 
This corresponds to lowering the slope of the theoretical line of Fig. SI This might ex- 
plain why the d In(w^) / dlnT measured at short times where we expect vibrations to be 
dominant lie on the lower side of the line. 

VI. DISCUSSION AND CONCLUSION 

The liquids studied show no universal glass-transition Lindemann criterion when we 
compare the MSD on the nanosecond time scale. Three of the liquids were studied at two 
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different pressures. They obey a pressure-dependent Lindemann criterion, as predicted 
by the elastic models. Thus the use of pressure reveals a connection which is probably 
masked by details in the molecular interactions and geometries when comparing dif- 
ferent liquids; this suggests the existence of an intrinsic Lindemann criterion for each 
substance at their pressure dependent glass transition. More extensive studies including 
several pressures along the glass-transition line of each liquid are needed in order to es- 
tablish the range in which this result holds. A failure of a universal Lindemann criterion 
might also be due to the fact that relaxations contribute to the MSD as was discussed in 
the previous section. Similar to the elastic model the Lindemann criterion is based on a 
vibrational picture. 

Above it appears that the elastic-model prediction underestimates the activation 
energy temperature dependence. We suggest that these deviations are caused by an- 
harmonic effects. In this context it should be noted that despite the simple "harmonic" 
appearance of the elastic models, anharmonicities must play a role, even in the simplest 
elastic models. Thus the effective, temperature-dependent elastic constant (or curvature 
at energy minima) reflects anharmonicity because in truly harmonic potentials the elastic 
constants are temperature independent. 

While the vibrational part of the mean-square displacement does not follow the predic- 
tion of the elastic models, we find that the total MSD measured, {u^){T), at the nanosec- 
ond time scale (vibrations and relaxations) approximately follows a proportionality of 
the type AE(T) a T / (w^), where AE(T) is the activation energy governing the alpha re- 
laxation. This one-to-one finding is based on measurements of MSD on the nano second 
timescale by neutron scattering as function of temperature of molecular liquids covering 
a significant range of fragilities. The proportionality AE(T) <x T / (u^) shows that there 
is a connection between the fast and the slow dynamics close to the glass transition. It is 
not clear how causal the relation is, whether the increase in MSD leads to higher mobility 
and consequently a speed up of the alpha relaxation, or the increased MSD is a due to 
a precursor of alpha relaxation itself, for example as a high frequency von Schweidler 
regime. 

To summarize we find (i) an intrinsic Lindemann criterion for each liquid as predicted 
by the elastic models by studying the same liquids at different pressures; (ii) the temper- 
ature dependence of the MSD on the nanosecond time scale links to the liquid fragility 
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as predicted by the elastic models. These observations are rationalized by introducing an 
anharmonicity in the elastic models, and they are fully consistent with other experimen- 
tal features and correlations found in the literature. The findings in this work underline 
that a full understanding of the viscous slowing down must involve both the fast and the 
slow dynamics, and suggest that elastic models offer a starting point for understanding 
the connection between these different time scales. 
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